Colloidal interactions in two dimensional nematics 
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The interaction between two disks immersed in a 2D nematic is investigated (i) analitically using 
the tensor order parameter formalism for the nematic configuration around isolated disks and (ii) 
numerically using finite element methods with adaptive meshing to minimize the corresponding 
Landau-de Gennes free energy. For strong homeotropic anchoring, each disk generates a pair of 
defects with one-half topological charge responsible for the 2D quadrupolar interaction between 
the disks at large distances. At short distance, the position of the defects may change, leading to 
unexpected complex interactions with the quadrupolar repulsive interactions becoming attractive. 
This short range attraction in all directions is still anisotropic. As the distance between the disks 
decreases their preferred relative orientation with respect to the far-field nematic director changes 
from oblique to perpendicular. 



PACS numbers: 77.84.Nh Liquids, emulsions, and suspensions; liquid crystals 
liquid crystal structure - 61.30.Jf Defects in liquid crystals. 
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I. INTRODUCTION 

In the last few years, there has been considerable theo- 
retical and experimental interest in colloidal systems with 
anisotropic host fluids Q. 

Since the temperature changes the anisotropic prop- 
erties of the host fluid, the dispersed particles exhibit a 
rich variety of collective behaviors, leading to unexpected 
phase diagrams B. Depending on the details of the sys- 
tem, chain- like [|3j,|4| , crystal || or cellular structures 
have been observed. 

In the case of nematic hosts in three dimensions (3D), 
the colloidal interactions result from the competition be- 
tween the anchoring properties of the spherical colloidal 
surfaces favoring e.g., radial (homeotropic) nematic ori- 
entation (anchoring) and the bulk elasticity that favors 
a uniform nematic. For sufficiently strong homeotropic 
anchoring W , the spherical colloids may induce topolog- 
ical defects or singularities in the nematic director field. 
If the size of the colloidal particle a is large when com- 
pared to the nematic correlation length £, a hyperbolic 
point defect appears at a certain distance from the parti- 
cle (satellite dipolar configuration Rflptfiol ) . For smaller 
particles, however, a ring disclination around the equato- 
rial plane of the sphere may have a lower energy (saturn- 
ring quadrupolar configuration pjj). Finally, if Wa is 
smaller than the typical nematic elastic constant K, the 
director field varies smoothly and there are no defects 
(surface-ring quadrupolar configuration Depend- 
ing on the symmetry of the distortio n |13[ , the dominant 
long range interactions are dipolar [|14| or quadrupolar 
[ plPl . The short range interactions are more difficult 
to predict, as nonlinear terms in the elastic free energy 
come into play. Experimental evidence from a system 



where silicon oil droplets were dispersed in a nematic 
host |l7j reveals that when the nematic configuration 
around each droplet has dipolar symmetry, the attractive 
long range interaction between aligned droplets becomes 
strongly repulsive at short distances. A nematic topo- 
logical defect appears between the particles and prevents 
their collapse. By contrast, when the long range dis- 
tortion exhibits quadrupolar symmetry no such defects 
appear and the oil droplets coalesce. 

2D colloidal systems were also investigated recently 
and similar qualitative features were observed. In these 
systems the interactions between the dispersed parti- 
cles are also mediated by the elastic distorsions of the 
anisotropic host. However, for 2D nematic hosts and 
homeotropic anchoring conditions, only quadrupolar con- 
figurations are stable |lq| . Dipolar configurations may be 
found in 2D smectic C films , a system that is similar 
to a 2D nematic, but where the director does not exhibit 
mirror symmetry, thus excluding configurations with half 
integer topological charges. In this case, Pettey et al. 
p(3f obtained analytic expressions for the free energy of 
a system of multiple disk-defect pairs, valid in the linear 
regime where the particles are sufficiently far apart. We 
investigated the short range interactions in this system 
by using a numerical method capable of describing the 
nematic orientation profiles induced by a small number 
of colloids with arbitrary shape in 2D |2l]] . Numerical in- 
vestigations of this type are faced with a problem due to 
the very different length scales characterizing the colloids 
and the topological defects. We have solved it pj[ , by 
using finite elements with adaptive meshing to minimize 
the Frank free energy and obtained very accurate inter- 
action energies between disks at arbitrary separations. 
We have found a strong repulsion at short range followed 
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by the expected long range dipolar interaction with a 
pronounced minimum at an intermediate distance. The 
equilibrium director field configuration exhibits a topo- 
logical defect at the mid-point between the disks. 

In this article, we use similar numerical methods pl| ] 
to minimize the Landau-de Gennes free energy of two 
circular disks in a 2D ncmatic. The great advantage of 
this model is that no special treatment of the topological 
defects is required, by contrast with the Frank elastic free 
energy description, where the defects appear as singular- 
ities of the director field n. In the next section, we de- 
scribe the nematic free energy based on the tensor order 
parameter formalism. Then, we generalize the analytical 
results of Pettey et al. |^| to nematic configurations with 
half integer topological defects and obtain analytic ex- 
pressions for the free energy of the quadrupolar nematic 
distortion for a single disk and the long range quadrupo- 
lar interaction between disks. Finally, we present the 
results of extensive numerical calculations for the inter- 
action between disks at arbitrary separations. 

II. THE LANDAU-DE GENNES FREE ENERGY 

In the tensor order parameter formalism the 2D ne- 
matic order is represented by the (tilt) angle of the di- 
rector field n = (cos 8, sin 8), that describes the mean 
orientation of the molecules, and the orientational order 
parameter Q, that measures the degree of alignment of 
the molecules with respect to the mean orientation. In 
the isotropic phase, where the molecules point in every 
direction, Q = 0. By contrast, in the nematic phase the 
molecules are preferentially aligned and the orientational 
order parameter takes a well defined equilibrium value 
Q = Qeq that depends on the temperature. 

We define a 2D symmetric traceless tensor field, from 
the tilt angle 8(r) and the orientational order parameter 
Q(r), 

Q ij (r)=Q(r)(n i (r)n j (r)-S ij /2) (1) 

that is invariant under the nematic symmetry operations 
n — > — n and r — > r. 

The Landau-de Gennes free energy may be understood 
as a Taylor expansion in terms of the tensor order param- 
eter Qij (the bulk free energy) and its derivatives dkQij 
(the elastic free energy) p2fl : 

F = L ( _ ^ + ? (Qlj)2 + k {dkQ ^ )2 Y r (2) 

where the summation over repeated indices is implied. 
The integration is over the 2D space occupied by the ne- 
matic fl. Only terms with the nematic symmetry are 
allowed in the free energy expansion. In addition, sta- 
bility requires that the total free energy is bounded from 
below. In the nematic A and C are positive constants, 



defining the equilibrium value of the orientational order 
parameter Q eq = ^JIAjC . For simplicity, we adopted 
the one-constant approximation for the elastic free en- 
ergy. In the general 2D case two elastic terms - corre- 
sponding to different invariant derivatives of the tensor 
order parameter - have to be included in the expansion. 
The elastic constant L is related to the Frank elastic con- 
stant through K = 4LA/C. The total free energy can be 
written as a functional of the tilt angle 8 and the orien- 
tational order parameter Q: 

F= L(~^ Q2+ S° 4 + i (VQ)2 + L Q 2{ y e ?Y*- 

(3) 

Let a be a typical length of the system, e.g. the radius 
of the colloidal particles. If L/a 2 <C A, the energy associ- 
ated with the variation of the tilt angle 8 is much smaller 
than the energy associated with the orientational order 
parameter. In this case we may set Q = Q eq everywhere 
in the nematic and the free energy simplifies to 

Fd= f L {ve)2d2r (4) 

which is the one-constant elastic Frank free energy. 

For strong distortions, however, topological defects 
may appear. In this case the director field becomes sin- 
gular and within a small region around the defect the 
orientational order parameter vanishes and the contribu- 
tion of the first three terms in Eq. (^|) has to be taken into 
account. This energy - the core energy of the defect - 
arises from a small region of size £, and may be estimated 
by a boundary layer analysis |23|| , 

Let us consider a topological defect of charge q placed 
at the center of the 2D space. In polar coordinates 
r = (r, (f>), the defect may be described by 8{cj)) = q<f>, 
asymptotically valid in the core region r < £, far from 
any other boundaries. The orientational order parame- 
ter, constant outside the core region, decreases rapidly to 
Q(r = 0) = inside the core in order to regularize the ne- 
matic free energy. A simple description of this behavior 
of Q is given by the ansatz 

Q(r) = Q eg (l - e~i) . (5) 

By substituting this expression into the free energy (|^) 
we obtain the simple result 

F(0 = F c (0 + F el (0 = f^ 2 + q 2 *K In ^ (6) 

where / w 0.97 and is the size of the system. Dif- 
ferentiation with respect to £ yields £ = |q|£, where 
Q = \J (2tt/ f)(L/ A) is the nematic correlation length. 
If C is much smaller than the other characteristic lengths 
of the system, we can use the Frank elastic free energy 
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(Q) to obtain the director configuration in systems with 
topological defects. The Frank elastic free energy, how- 
ever, is defined outside the core region only. The total 
free energy can be computed by additing to it the core 
energy of each defect 



,TtK 



(7) 



III. ISOLATED DISK 

Fukuda and Yokoyama [^8) studied the distortion of a 
2D nematic around an isolated disk, with homeotropic 
anchoring, by minimizing the Landau-de Gennes free en- 
ergy (^) using numerical methods with adaptive grids. 
They have shown that a director configuration with a 
dipolar defect (topological charge q = —1) is unstable, 
and changes into a configuration with a pair of defects 
with topological charge q = —1/2 that are located sym- 
metrically around the particle. 

The energy of this quadrupolar configuration may be 
calculated analitically and compared with the energy of 
the dipolar case, obtained by Pettey et. al. |Q for 2D 
smectic-C films. 
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FIG. 1. Single disk. Left: notation and reference frame 
used in this work. Right: quadrupolar configuration with two 
half integer defects (grey) . The lines represent the orientation 
of the nematic director. 

Consider the disk at the center of the reference frame 
r = (x, y), or r = (r,<fi) in polar coordinates, where <j> 
is defined counterclockwise, starting from the x-axis, in 
the range — 7r < <fi < ir (see Fig{l]). At the outer bound- 
ary, for large r, the nematic is uniform with 9 = tt/2. 
Homeotropic boundary conditions apply at the disk's 
boundary, i.e., 0(<f>) — (f> on the upper contour and 
0((f>) — 7r + on the lower one. Minimization of the Frank 
elastic free energy ([I]) yields Laplace's equation AO = 0, 
everywhere (in the nematic) except within the defects. 

An image solution with two half integer topological 
charges located symmetrically at a distance from the 
center of the disk, satisfying the homeotropic boundary 
conditions at the disk's contour and a uniform director 
at infinity, is given by 



= \ - 2^ + \A d + \i>-r d + \i>f d + \lp-f d (8) 

where = a 2 jr^ is the position of a virtual half integer 
topological defect inside the disk. The field 
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describes a topological defect of charge q = 1 located at 
r i = (ccj , 0) and satisfies Laplace's equation everywhere 
except at r^. 

The elastic Frank free energy can be calculated using 
Green's Theorem: 



I (V0) 2 d 2 v = I 9^-dl - / 9A9d 2 r , 
Jn Jan dn Jn 



(10) 



where 90 is the boundary of the 2D nematic domain, 
dl = {dx, dy) is an oriented contour element, and 
d9/dn = {-d9/dy,d9/dx). If the tilt angle 9 satisfies 
Laplace's equation, the second term on the right hand 
side of vanishes, and the free energy becomes 
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Close to the outer boundary, for r 3> a, the director 
is nearly uniform and the contribution of the integral 
over this contour may be neglected. The contribution of 
the integral over the upper disk's contour is calculated by 
taking 7i(i) = (x(t),y(t)) — (acost, asmt) with 6(t) = t, 
yielding 
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By symmetry the integral over the lower disk's contour 
is identical. 

Now consider the integrals over the straight lines along 
the x-axis from the disk to the defects. On 72 and 76, 
the tilt angle is zero and there is no contribution to the 
free energy. On 74 (t) = (t, 0) (and 7s), 9(t) = tt and a 
simple calculation yields 
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(13) 



in the limit rd >> £. If a » £ we may neglect the con- 
tributions from the circular contours around the defects. 
Finally, the Frank elastic free energy is given by 
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Minimizing with respect to r d yields (r<j) eg = y/7/3 a, 
as obtained by Fukuda and Yokoyama using a sim- 
ple force balance argument. The total free energy for 
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this quadrupolar configuration is obtained by adding the 
core energy of the defects 2F C — irK/4, and using the 
appropriate core's size £ = (/2. 

The Frank elastic free energy for a dipolar configura- 
tion, with a topological defect of charge q = — 1, with 
the preferred position at = (0, Td) was calculated by 
Pettey et. al. pof: 
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eq = V2 a. The 
= nK/2, and the 



yielding the equilibrium position (r<j 
core energy of the defect is now F c 
core's size is £ = (. 

The difference between the free energy of these config- 
urations, when the core energies of the defects are taken 
into account, is a function of the correlation length 



Fdip — F qua d ~ — — ^2.08 + In — ^ 



(16) 



For a > C, the difference is positive, and thus the 
quadrupolar configuration has the lowest free energy. 

This result contrasts with that for spherical colloidal 
particles in 3D nematic liquid crystals. In the latter case, 
the dipolar configuration has a punctual defect (with zero 
core energy). The quadrupolar configuration exhibits a 
disclination ring, and its elastic free energy depends log- 
arithmically on the ratio between the correlation length 
and the particle's radius. For small particles, the nematic 
distortion is of the quadrupolar type. As we increase 
the ratio a/£, the energy of the quadrupolar configura- 
tion increases, and eventually exceeds that of the dipolar 
configuration. This crossover occurs at et/£ ~ 10 3 , that 
corresponds to particles with radii of the order of /im. 



IV. LONG RANGE INTERACTION 



If the separation R = (i?cosa, i?sina) between parti- 
cles is large, the nematic distortion will be approximately 
the sum of the isolated quadrupolar solutions, w 61+62, 
with an interaction energy 
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where the subscripts refer to the different disks (see 
Fig.2). 




FIG. 2. Two particle configuration: notation and reference 
frame used in the analysis. 

Consider the integral over the contour T2 around the 
second disk. We use the same notation as for isolated 
disks and neglect the contributions from the circular con- 
tours around the defects. To calculate the contribution 
from the upper circular contour around disk 2 we write 
71 (t) = (i? cos a + a cos t , R sin a + a sin t) . On this con- 
tour 62 (t) = t, and 6\ can be expanded as a power series 
in a/R yielding 
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When the contribution from the lower circular contour 
around the disk is added, some of the terms in ([l8|) can- 
cel. After integration, the sum of these two contributions 
gives 
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The integrals over the straight lines can be computed in 
a similar manner. On 72 and 76 the tilt angle 82 = 0, 
and there is no contribution to the free energy. On 
74,{t) = (Rcosa + t, Rs'ma) (or -f S ), 6i{t) — ir and Q\ 
may be expanded as a power series in a/R. The result is 
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The derivative d 2 81 /dxdy is obtained from the 
quadrupolar expression 
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valid for large r = y 1 x 2 + y 2 . Finally, adding the sym- 
metric contribution from the contour Ti, the total inter- 
action energy becomes 
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The effective interaction between disks decays as i?~ 4 
and is strongly anisotropic: repulsive if the particles are 
aligned horizontally or vertically (a = or a — tt/2, re- 
spectively), and attractive for intermediate oblique ori- 
entations (the preferred orientation is a — 7r/4). 



V. SHORT RANGE INTERACTION 



At small disk separations, the nematic deformation is 
no longer the sum of the isolated quadrupolar solutions 

The interaction between two circular disks was studied 
numerically, using finite elements with adaptive meshes 
to minimize the total nematic free energy at a given sep- 
aration. A first triangulation respecting the predefined 
physical geometry (disks separation and orientation) is 
constructed j24j] . The triangulation space corresponds to 
the domain of integration of the free energy. It is de- 
limited by the disks' boundaries and an outer boundary 
at infinity, in practice a large distance from the disks, 
of the order of 20a. Although analytically simpler, the 
Frank free energy (^) is plagued with divergences of V0 
at the defects, and thus it is not the most adequate for 
numerical computations. For this reason, we used the 
Landau-de Gennes free energy (||) for a 2D nematic. The 
functions Q xx and Q xy are set at all vertices of the mesh 
and are linearly interpolated within each triangle. Us- 
ing standard numerical procedures the elastic free en- 
ergy is minimized under the constraints imposed by the 
boundary conditions, i.e., homeotropic strong anchoring 
at the disks perimeter and constant alignment at the 
outer boundary. Finally, to overcome the differences in 
length scales set by the disks and the defects, of up to 2 
orders of magnitude (we chose £ 2 ~ L/A = 10~ 3 ), new 
adapted meshes are generated iteratively from the result 
of the previous minimization. The new local triangle sizes 
are related with the variations of the previous solution. If 
the variations are strong, finer meshes are required in or- 
der to guarantee an almost constant numerical weight of 
each minimization variable The final meshes with 

2 x 10 4 variables and a minimal edge length of 10~ 4 a 
provided excellent accuracy, yielding free energies with a 
relative error of the order of 1CP 4 . 

To test the numerical accuracy of the procedure we 
calculated the nematic configuration for a single isolated 
disk, and found a solution with two half integer defects 
at rd = 1.23a ± 0.01a, in very good agreement with the 
analytical result. 
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FIG. 3. Reduced elastic free energy F — F/K as a func- 
tion of the distance R between the disks, at three different 
orientations (a = 0, tt/4, tt/2). F u = F[Q = Q eq \. 

The effective interaction between two disks is plotted 
in Fig. 3 as a function of the distance R between the disks, 
at three different orientations (a = 0, 7r/4, 7r/2). A sim- 
ple fit of these curves confirmed the quadrupolar long 
range decay R , at large disk separations. However, 
for smaller separations, the free energy changes dramat- 
ically. In particular, orientations that are repulsive at 
large separations (a = 0,7r/2) become attractive at dif- 
ferent threshold distances (Rth ~ 5a and Rth ~ 3a, re- 
spectively). Also, the free energy corresponding to the 
parallel orientation (a = 0) becomes smaller than the 
free energy of the long range preferred oblique orienta- 
tion (a = 7r/4) at short distances. 

To elucidate this behavior we plotted on the left 
of Fig. 4 the elastic free energy as a function of the 
orientation a, at three different separations (R — 
2.4a, 3.0a, 4.0a). The minimum of each curve defines a 
preferred orientation a*, that decreases as the disks ap- 
proach each other. On the right of Fig. 4 we see that at 
very large separations, the disks have an oblique orienta- 
tion (with a* = tt/4), that changes to a parallel orienta- 
tion (with a* — 0) as their separation decreases. 
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FIG. 4. Left: Reduced elastic free energy as a func- 
tion of the orientation a, at three different separations 
(R = 2.4a, 3.0a, 4.0a). Right: Preferred orientation a* as 
a function of the distance R. 



In the following figures we show the nematic order 
parameter for several separations, when the orientation 
of the disks is parallel (a = 0, Fig. 5) or perpendicular 
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(a = 7r/2, Fig. 6). The nematic order parameter varies 
between Q = (grey regions) and Q = Q eq (white). 
These regions around the disks correspond to the nematic 
defects. 




FIG. 5. Nematic configurations for several separations and 
horizontal alignment a = 0. The nematic order parameter 
varies between Q = (grey regions) and Q — Q eq (white). 



lies between 30 < W/Ka < 40. In the weak anchoring 
regime the coalescence of the dropplets may occur. 




a R/a b " R/a 



FIG. 7. Reduced elastic free energy as a function of the 
disks separation R, for three reduced anchoring strengths 
(W/Ka — 1000; 40; 30), corresponding to the circles, dia- 
monds and triangles respectively. Left: Parallel alignment 
(a — 0). Right: Perpendicular alignment (a — 7r/2). 



In both disk alignments, if the separation is large 
enough, the half integer defects around each disk stay 
at the same positions as in the isolated disk configura- 
tion. However, at a certain separation Rth, exactly when 
the interaction becomes attractive for that orientation, 
the position of the defects starts to change, in order to 
minimize the total free energy. 




FIG. 6. Nematic configurations for several separations and 
perpendicular alignment a — n/2. 

Finally, we analyzed in detail the interaction between 
the disks for very small separations. For strong (fixed) 
homeotropic anchoring, we observed a repulsion when 
the disks are nearly at contact (for all orientations) at a 
separation R ss 2.1a. This repulsion may not occur if we 
do not impose fixed homeotropic boundary conditions at 
the disks. If we consider a finite anchoring strength we 
have to add to the free energy (|J) a surface term 

F a = \ I (Qij - Q^fdl (23) 

Z JAi+A 2 

where Q™ are the values of the tensor order parameter 
for homeotropic aligncmcnt. 

In Fig. 7 we plot the total elastic free energy as a func- 
tion of the distance R between disks, for three differ- 
ent reduced anchoring strengths (W/Ka — 1000; 40; 30). 
These anchoring strengths are large enough to induce 
defects in the nematic host. Clearly the repulsion ap- 
pears at small separations for strong anchoring (W/Ka = 
1000) and vanishes at a critical anchoring strength that 



VI. CONCLUSIONS 

We compared the (analytic) elastic free energies of 2D 
nematic configurations with one (satellite) and two (ring) 
defects around isolated disks with homeotropic boundary 
conditions, and found that the ring configuration is the 
stable one for realistic values of the reduced correlation 
length (/a. 

Using the analytic ring solution we derived the long 
range interaction between two disks immersed in 2D ne- 
matic hosts. We found an anisotropic interaction, that is 
repulsive for orientations close to perpendicular or par- 
allel orientations and attractive for oblique aligncmcnt, 
with a preferred orientation at a = 7r/4. In this limit, 
the interaction decays as i? -4 . 

Finally, the analytic results were compared with very 
accurate numerical calculations, valid at arbitrary sepa- 
ration. At small disk separations, the defects may change 
their positions, leading to an unexpected complex inter- 
action with the long range repulsive orientations becom- 
ing attractive. As a result as we decrease the distance be- 
tween the disks, their preferred orientation changes from 
oblique to parallel alignment. For sufficiently large an- 
choring strengths, at very short range, we found a new re- 
pulsion for all orientations, that may prevent coalescence. 
This repulsion vanishes for small anchoring strengths. 

The equilibrium structure of inverted nematic emul- 
sions is complex and depends sensitively on the radii of 
the disks, the nematic correlation length and anchoring 
strength. In particular short range nonlinear effects may 
change the interaction between disks dramatically. It 
is thus reasonable to assume that three body interac- 
tions may be relevant to describe the formation of the 
the structures observed experimentally. Finally, in order 
to understand the observed structures dynamical aspects 
may also have to be taken into account. 
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